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Electrons, when scattered by static random disorder, form standing waves that can be imaged 
using scanning tunneling microscopy. Such interference patterns, observable by the recently de- 
veloped technique of Fourier transform scanning tunneling spectroscopy (FT-STS), are shown to 
carry unique fingerprints characteristic of the electronic order present in a material. We exploit this 
feature of the FT-STS technique to propose a test for the nature of the enigmatic pseudogap phase 
in the high- Tc cuprate superconductors. Through their sensitivity to the quasiparticle spectra and 
coherence factors, the FT-STS patterns in principle carry enough information to unambiguously 
determine the nature of the condensate responsible for the pseudogap phenomenon. In practice, 
the absence of a detailed understanding of the scattering mechanism, together with the experimen- 
tal uncertainties, prevent such an unambiguous determination. We argue, however, that the next 
generation of FT-STS experiments, currently underway, should be able to distinguish between the 
pseudogap dominated by the remnants of superconducting order from the pseudogap dominated by 
some competing order in the particle-hole channel. Using general arguments and detailed numerical 
calculations, we point to certain fundamental differences between the two scenarios and discuss the 
prospects for future experiments. 



I. INTRODUCTION 

The problem of high- Tc superconductivity continues to 
engage the scientific community, yet much remains un- 
known. The phase diagram of the cuprates is being elu- 
cidated with unprecedented accuracy but a unified mi- 
croscopic theory of the observed phases is still lacking. 
While the question of how exactly the Mott-Hubbard 
antifcrromagnetic insulator evolves into a high- Tc super- 
conductor upon hole doping remains unanswered, a num- 
ber of phenomenological approaches to the problem have 
been developed. In the absence of a microscopic theory, 
progress can be made by means of a phenomenological 
description of the various phases, some of which are mo- 
tivated by microscopic argumentsiiiSiiii^ 

Among the various phases, two are especially well es- 
tablished: the antiferromagnet (AF) near the half fill- 
ing and the superconducting phase which occurs upon 
hole doping. The latter is reasonably well described by a 
BCS-like theory with a d-wave symmetric order param- 
eter, Ak- It is this symmetry that gives rise to unique 
physics of the low energy excitations which reside in the 
vicinity of the nodes of the gap function. 

Intermediate between the AF and d-wave superconduc- 
tor (dSC) is the so called pseudogap phase which exhibits 
a gap in the single particle density of states but is non- 
superconducting.'' Experimentally, the pseudogap phase 
is characterized by a d-wave gap in the electron excita- 
tion spectrum,- which sets in at a temperature T* much 
higher then the superconducting transition temperature 
Tc. The low energy nodal quasiparticles seem to exist in 
the pseudogap phase as well, however, they are strongly 
interacting and exhibit a high scattering ratei^ 

Various authors have made proposals to describe the 
pseudogap phase using different phenomenological ap- 



proaches. One can identify two broad classes of theo- 
ries. One class includes theories of competing orders, in 
which the pseudogap results from the formation of static 
or fluctuating order in the particle-hole channel. Exam- 
ples include the AF spin density wave (SDW) , the charge 
density wave (CDW), stripes, and the so-called d-density 
wave (DDW, also known as the flux phase). This order 
parameter is assumed to compete with superconductivity 
and its formation is energetically favorable in the pseudo- 
gap phaseASiifliii The other class of theories follow the 
ideas of order parameter phase fluctuations, first intro- 
duced by Emery and Kivelsoni This school of thought 
views the pseudogap as a superconductor on short length 
scales that fails to superconduct on longer length scales 
due to phase fluctuations. In particular, this implies that 
a local superconducting order parameter is formed at T* , 
but phase coherence is only achieved at y^jiiSiiSiiSiiiiiii^ 
Unraveling the nature of this enigmatic phase is one of 
the central issues in the physics of high- Tc cuprates. 

In this Paper we elaborate on a test for the nature of 
the pseudogap phase using the Fourier transform scan- 
ning tunneling spectroscopy (FT-STS), proposed by us 
earlier We demonstrate that the new generation of FT- 
STS experiments, performed in both the dSC and pseu- 
dogap phases, can provide a clear distinction between 
the two types of order, particle-hole (p-h) or particle- 
particle (p-p), and can therefore serve as a referee be- 
tween the above two classes of theories. Such a distinc- 
tion is possible due to fundamental differences between 
the two types of order. In a superconductor the elemen- 
tary excitations consist of an admixture of electrons and 
holes. The charge of these quasiparticles is not a good 
quantum number, and each quasiparticle carries a quan- 
tum phase that describes its particle-hole mixing. As a 
result, scattering processes which give rise to the FT- 
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STS signal, have high or low probability depending on 
the relative phase of the scattering quasiparticles. This 
probability is given by the relevant coherence factors. On 
the other hand, when the scattering amplitude may de- 
pend on various microscopic details, we show that there 
exist certain generic features which only depend on the 
type of ordering present in the system. 

The most straightforward test follows from the fol- 
lowing simple consideration. The interference patterns 
caused by a p-p order parameter are known; these are 
just the experimental results observed in cuprates below 
Tc^i^ Comparing these patterns with the experimental 
results in the pseudogap phase should reveal the nature 
of the pseudogap order parameter. If the patterns are 
similar to those of the dSC state, it is likely that the 
pseudogap order is in the p-p channel. If, on the other 
hand, an abrupt change in the patterns is found, the con- 
clusion will be that a different type of order, in the p-h 
channel, governs the pseudogap. We shall discuss the 
relation to the available experimental data in Sec. IV. C. 

This article is organized as follows. In the next section 
we revisit the formalism of the quasiparticle interference 
in the superconducting state, reviewing some recent the- 
oretical work and commenting on the agreement with 
experiment. We choose a framework for further analy- 
sis and discuss how it may be justified on microscopic 
grounds. In section III we apply the same theoretical 
treatment to the pseudogap phase. We give a general 
discussion of the problem and then study two represen- 
tative theories of the pseudogap, the QED3 theory and 
d-density wave (DDW). We comment on their relevance 
in light of experimental observations. Our conclusions 
are given in section IV. 



II. THE SUPERCONDUCTING STATE 

A. General considerations 

In the high resolution FT-STS experiments the spec- 
trometer's sharp tip is placed above a freshly cleaved 
sample surface, and a bias voltage is applied between the 
sample and tip. Tunneling occurs through the vacuum 
barrier and produces electric current which is measured 
as a function of the bias voltage. The measured differen- 
tial conductance dI{V, r)/dV is proportional to the local 
density of states (LDOS) n(r, oj) of the sample below the 
STS tip at the point r. A wide field of view of about 
600Ax600A is scanned with atomic resolution, and at 
each bias voltage the spatial map of n{r,uj) is collected. 
The spatial Fourier transform of this quantity, n(k, uj), is 
then studied. The resulting fc-space distributioni2ii2iiS 
reveals distinctive patterns with peaks at special wave 
vectors that correspond to certain periodic structures in 
the real space electron wavefunctions. In a clean, homo- 
geneous superconductor such structures would be absent; 
they apparently result from the translational symmetry 
breaking caused by disorder. At present, the details of 
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FIG. 1; Top right: The banana shape contours of constant 
low energy around the nodes of the gap and seven vectors 
connecting their ends at energy = 0.7Ao^ Top left and 
bottom: Numerical evaluation of the Fourier transformed lo- 
cal density of states in Eq. H12|l with t' = 0.3t, Ao = 0.2t, 
fj, — —t and u = 0.35, 0.5 and 0.7Ao respectively. In order 
to best fit the experimental datai^*— we present the real part 
of A(g, oj) and suppress the A+A- term as explained in the 
text. 



this scattering potential are not fully understood. 

These pioneering FT-STS experiments were followed 
by several theoretical studies that calculated the LDOS 
modulations in the BCS c?-wave model using various 



approximations 
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First, a heuristic picture re- 



ferred to as the octet model was put forward by Hoff- 
man et a/ii2*i242& and provided a simple and appealing 
explanation of the data. The octet model identifies eight 
points in momentum space at which there is the largest 
single particle density of states (DOS), and asserts that 
the experimental FT LDOS is expected to have a peak 
at any point q in the Brillouin zone such that q connects 
any two of the octet points. As illustrated in Fig. ^ the 
latter occurs at the ends of the banana shaped contours 
of constant energy (CCE), where the energy is equal to 
the bias voltage, hui = eV. 

The original octet model explains the experimental 
data surprisingly well, considering that it ignores several 
important issues. In particular, the on-shell excitations 
away from the octet points are not taken into account. 
One can argue that such excitations have lower single- 
particle DOS than the octet model points. However, the 
ratio between the single particle DOS at the edges of the 
CCE bananas to the DOS of any other point on the CCE 
is only a result of the velocity anisotropy and is not as 
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significant as other considerations. This can be seen by 
expanding the energy dispersion near the nodes and cal- 
culating the DOS. In the "nodal" approximation, scaling 
the momentum coordinates by the appropriate velocities 
{vp and va) leads to a circular CCE, on which there is 
no special point that could define a peak position. 

A more detailed analysis compares the measured spec- 
trum to a joint density of states (JDOS), D{q,uj), given 
by a k-space convolution of the single particle density of 
states 

i?(q,c^) =^5(c^-Sk)'5(w-^k-q), (1) 
fc 

where is the single particle excitation energy in the 
superconductor. This approach includes the on-shell pro- 
cesses away from the nodal points but neglects the ofF- 
shell processes as well as the coherence factors. We shall 
see shortly that in a d-wave superconductor the coher- 
ence factors for scalar disorder are precisely such as to 
greatly enhance the contributions from the octet points. 
We also find that the off-shell processes mostly contribute 
to the background and, interestingly, tend to worsen the 
agreement with experimental data. 

It is interesting to note that the theoretically more sat- 
isfying approachesiSi22i2ii22., which model the LDOS us- 
ing the full electron Green's function (computed within 
various reasonable approximations), generally fail to ac- 
count for the details of the experimentally observed pat- 
terns. The main problem is that these theoretical ap- 
proaches typically produce a number of prominent back- 
ground features that are not present in experiment. As 
a result, the computed patterns do not really resemble 
the experimental data (compare, e.g., Refs. I20I and ll^ . 
even if they contain many of their characteristic features. 
At present it is not clear whether this discrepancy is a 
sign of some fundamental difhculty in the existing mod- 
els or simply a question of details that have not been yet 
properly addressed. 

Here we take the latter view; i.e. that a simple BCS- 
type model of a d-wave superconductor with appropriate 
disorder should provide a complete description of the FT- 
STS data once the details of the scattering mechanism 
have been ironed out. In the present paper we are inter- 
ested in these details only to the extent that will allow 
us to establish a baseline for our discussion of FT-STS 
in the pseudogap state. With this goal in mind in what 
follows we design a model for scattering in the dSG that 
produces patterns in reasonable agreement with the ex- 
periment. We also explore several routes to justify this 
model from microscopic considerations, but in this we are 
only partially successful. We then take this same model 
and use it to discuss the pseudogap state. 

B. The model 

We shall describe the system by a standard BCS d- 
wave Hamiltonian with disorder, treated as a perturba- 



tion. In the following, we consider three types of dis- 
order, namely in the charge, spin and pairing channels. 
We show how different types of disorder produce differ- 
ent scattering patterns and comment on their relation 
to the patterns observed in the FT-STS experiments. 
Our starting point is similar to that used by previous 
authors22iSi42^ but there will be some differences which 
will allow us to closely reproduce the experimental data. 

The system is described by the Hamiltonian H = Ho + 
SH, where 

is the d-wave BCS Hamiltonian describing the pure, uni- 
form dSC. 4'k ~ ("^kT' '^-ki) is the Nambu space creation 
operator and ek , Ak are the band structure and gap func- 
tion respectively. In what follows, we choose to work with 
the second-nearest neighbor hopping model and a dx^-y^ 
gap function 

Ek = — 2t{cos kx + COS ky) — 4,1' COS kx COS ky — fi, 
Ak = ^Ao{cos kx - cos ky). (2) 

Disorder is generically described by a Hamiltonian 

'^^= E = E^k^kk'^k'. (3) 

i.j,(T,(7' k,k' 

Here i,j is a lattice index, V^"^ is the random potential 

and Vkk' its Fourier transform, a 2 x 2 matrix in Nambu 
space. 

The electron propagator of the perturbed system can 
be expressed as 

G(k, k',uj) = G°(k, c^)5k.k' + G°(k, c^)fkk'(c^)G°(k', uj), 

(4) 

where 

G°(k,c^) = [c.-a3ek-^iAk]-i (5) 

is the bare Green's function and Ui are the Pauli matri- 
ces. 'rkk'(w) is the T-matrix, subject to the Lippman- 
Schwinger equation 

Tkk'(w) = Vkk' + J2 VlcqG°(q, w)Tqk'(c^). (6) 
q 

The local density of states at the point r of the sample 
is then given by 

n{r,uj) = -ilm[Gii(r, r,w) -f G22(r, r, -w)], (7) 

where G(r, r',ij;) is the spatial Fourier transform of 
G(k, k',w). Here and hereafter "Im" is understood as 
taking the usual difference between the retarded and the 
advanced quantity measuring the discontinuity of the 
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propagator across the real frequency axis. It is impor- 
tant to keep this in mind since Vkk' is in general complex 
but its imaginary part does not contribute to the above 
discontinuity and therefore to LDOS in Eq. Q. Taking 
the "Im" symbol literally would lead to incorrect results. 

For arbitrary disorder one must typically solve for the 
T-matrix from Eq. (jSJ using numerical methods. We 
shall not attempt such a detailed numerical solution here. 
Rather, we shall limit ourselves to two different approxi- 
mations where progress can be made analytically. First, 
we consider the limit of weak but otherwise arbitrary dis- 
order. Second, we study the case of arbitrarily strong but 
dilute point-like scatterers. 



C. Weak scattering (Born limit) 




To leading order in the scattering potential V, i.e., in 
the Born approximation, we have T'kk'(w) — 14;k' and the 
perturbed Green's function is given by 

G(k,k',c^) = G°(k,c^),5k.k'+G°(k,c^)-l/kk'G"(k',c^). (8) 

Since the first term of Eq. (jSJ is uniform, all the inter- 
esting information resides in the second term. 

We now focus on the disorder potential in the charge 
channel. In this case we have on-site random potential 
Vi coupled to electron charge. 



(9) 



k,k' 



where Vq is the Fourier transform 



i.e., Vkk' = CTsVk-k 
of Vi. The advantage of the Born limit is that one may 
express the interesting non-uniform part of the measured 
FT LDOS, 5n(q, w), in a factorized for"^ 



Sn{q, uj) = --\Vq\lra[Aniq,u;) + A22iq, -Lu)] , (10) 



where 



A(q, u;)^Yl ^"C^' ^)<^3G"{k - q, cu) 

k 



(11) 



As discussed by Capriotti et al. Vq is a random func- 
tion with no interesting structure while A(q, lj) contains 
information about the underlying clean system. Any 
peaks observed in (5n(q, lu) therefore must be attributed 
to A(q, w). Thus, by studying the interference patterns 
in a disordered sample one can learn about the properties 
of the system in the absence of disorder. 

With the d-wave BCS Green's function, Eq. lO, one 
finds 

... X 1 (ia; + £+)(za; + £-)- A+A- 
Aii(q,-) = ;^L + + , (12) 



Al. 



with e± = ek±q/2, A± = Ak±q/2 and E± = ^ e 
Eq. H12|) can be evaluated exactly using numerical meth- 
ods. However, before doing so we shall elucidate its main 



FIG. 2: The nodes of the dSC gap on the underlying Fermi 
surface with the velocities vp and va as defined in each node. 
The vectors Q and Q' span antipodal and adjacent nodes 
respectively. 



features by employing the nodal approximation. In this 
approximation, valid at energies low compared to the 
maximum gap Aq, the spectrum is linearized in the vicin- 
ity of the nodes, and we consider scattering processes ei- 
ther within the same node or between different nodes. 
Let us first evaluate the LDOS resulting from intra-node 
scattering, which corresponds to the octet vector qyi^ 
For this case, Eq. H12|l can be approximated in the follow- 
ing way: near the node we use the 45° rotated coordinates 
ki and fc2 shown in Fig. |(2Jl. To first order in the momen- 
tum we expand ek ^ vpki and Ak — ^ WAfe- We then 
rescale the momentum coordinates so that vpki —^ fci, 
VAk2 k2, and redefine ^vpqi qi and ^i;a'?2 — > 92 • 



In the scaled coordinate frame 
simplifies to 



|k| and Eq. ^ 



A 



1 



lin 



d^k 



-OJ 



+ {kl-kl)^{ql-ql) 



VFVA J (27r)2 [a;2 + (k + q)2] [^2 ^. „ q)2i 



(13) 

As shown in appendix ^ the explicit solution can be 
found using the Feynman parameterization^^ and reads: 



Aiin(q,t^) 



1 



2itvfva 



x/i2~T 



: arctan ■ 



1 



(14) 



The function J-{z) implies an inverse square root singu- 
larity in both the real and imaginary parts of Aiin(q, a;) 
along an elliptic contour of constant energy given by 
Eq = 2lu (here Eq — 2q = \Jv\q{ -|- v\q?^). More impor- 
tantly, this singularity is weighted by an angular factor 
(92/9)'^ = {vA(l2lEqf, producing the largest amphtude 
at the two ends of the ellipse, as illustrated in Fig. |j3^). 
These points of largest intensity coincide with iqy. It 
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is important to note that the angular dependence results 
purely from the coherence factors (i.e. terms in the nu- 
merator of Eq. [Qjl and would not appear in the JDOS 
Eq. CJ. 

This particular angular dependence can be understood 
by analyzing the relevant BCS coherence factors as fol- 
lows. Consider first q that connects the two CCE point 
along the nodal direction. The process represents scat- 
tering of a pure particle to a pure hole since along this 
direction there is no p-h mixing due to the vanishing gap. 
The amplitude of such a process is zero due to the charge 
coupling. On the other hand the two points at the edges 
of the CCE both have equal particle and hole mixing 
(|uk| = kk| in the BCS notation) and the same quan- 
tum phase. This leads to constructive interference and 
maximal A(q, w). 

Similarly, we can treat magnetic disorder. In this case, 
the potential couples to the electron spin and is therefore 
proportional to the unit matrix ao=ll in Nambu space: 



5n 



mag r t 



k.k' 



(15) 

Again, one can describe the scattering patterns using 
A(q, w). It is easy to see that the expression for this 
quantity is just like Eq. H12() except for the sign change 
in front of the A_|_A_ term. In the nodal approximation 
for intra-nodal scattering we thus get 



vfVaJ (27r)2 [c^2 + (k + q)2][^2 + (k_q)2i 



The solution (see Appendix A) is now very different, 



(16) 



-1 



2'KVpVA 



J-'{z) — \J ~\ arctan 



1 



1 



(17) 



Here, since the spin is a good quantum number, no an- 
gular dependent coherence factors are present, and the 
interference amplitude is constant along the elliptic con- 
tour 2u} — \/vpqf + v\q^, as shown in Fig. (|3i). The 
square root singularity is replaced by a cusp along the 
contour. This result provides another example of the im- 
portance of the coherence factors. 

The linearized treatment for internodal magnetic and 
non-magnetic scattering is similar and given in ap- 
pendix A composite picture including scattering 
processes between all nodes is shown in panel (e) of 
Fig. Q. We observe that it indeed contains peaks and 
qualitatively exhibits features observed in FT-STS ex- 
periments. The linearized treatment also indicates that 
non-magnetic scattering alone cannot produce all of the 
experimentally observed features. The octet points qi, 
qs and q2,6j corresponding to adjacent nodes scattering, 
can all be produced together only by including magnetic 
scattering (see Fig. (|St)). This combination of magnetic 




FIG. 3: Interference patterns in the nodal approximation, a) 
Non- magnetic scattering within a single node Eq. b) Non- 
magnetic scattering between two antipodal nodes, Eq. ljA4^ . 
c) Magnetic scattering between adjacent nodes, Eq. 1A5II . d) 
Magnetic scattering within a single node, or antipodal nodes 
and non magnetic scattering of adjacent nodes. In (a,b,d) we 
used vf/va ~ 4 and in (c) vf/va = 1- Note that the axes 
here represent the nodal coordinates, qi (horizontal) and q2- 
e) A composite of the intra- and antipodal-node terms for 
non-magnetic scattering and adjacent-node magnetic scatter- 
ing showing the full Brillouin zone distribution in the nodal 
approximation. 



and non-magnetic scattering is discussed in the next sub- 
section. 

Beyond the linearized approximation, the interference 
patterns can be evaluated exactly, by numerical meth- 
ods, using the full model dispersion Eq. Work along 
these li nes has be en done for a single point like scatterer 
m Refs. 12(12 J2I and for arbitrary weak disorder in Refs. 
iUiSl . To evaluate Eq. ((T^ numerically we find it con- 
venient to perform the momentum space Green's func- 
tion convolutions by employing the fast Fourier transform 
(FFT) technique. In particular we transform G''(k,w) 
to real space, multiply the appropriate components, and 
perform an inverse FFT back to the momentum space. 
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FIG. 4: Top six panels: numerical evaluation of A{q,Lj) for 
ui = O.32A0, where the scattering potential is, from left to 
right proportional to 0-3 (non-magnetic scattering), ctq = 1 
(magnetic scattering) and an equal combination of both. Bot- 
tom: The octet points at the edges of the banana shaped 
CCE (left) and the experimental data of 14meV« 0.32Ao 
from Ref. (right). Note that the Brillouin zone of the 
experimental data is 45° rotated with respect to the crystal- 
lographic directions. 



Since the run time of FFT is only o(n log(7T.)) where n is 
the array size, this allows for fast evaluation of A(q, us) 
on large lattices (of up to 1024 x 1024 lattice points). 

The comparison of numerical results with the exper- 
imental data shown in Fig. Q leads to two interesting 
insights. When considering the Born limit of point-like 
charge/spin impurities represented by 173 or ctq matri- 
ces respectively, it appears that the real rather then the 
imaginary part of A(q, to) is describing the experiment 
more accurately. In addition, one achieves even bet- 
ter agreement with the data when neglecting the off- 
diagonal part of the unperturbed Green's function [i.e. 
A+A_ term in Eq. H12|) ]. Operationally, such suppres- 
sion of the off diagonal Green's function may be achieved 
as a result of an even combination of magnetic and non- 
magnetic scattering since the two contain contributions 
proportional to A+A_ with opposite signs. Fig. Q 
shows our numerical calculation of A(q, lu) for a mag- 
netic, non-magnetic and an even combination of both 
types of scattering. Our best fit to the data is obtained 
by taking the real part of the equal mixture of magnetic 
and non-magnetic scattering and is presented in Fig. Q . 



We adopt this result for FT-STS in dSC as our base- 
line for studies of the pseudogap state. Before we proceed 
to study the latter we briefly discuss how one could jus- 
tify this model based on the physics beyond the Born 
approximation. 



D. Beyond the Born limit: the T-matrix 

As mentioned above, the Born limit scattering model 
is in best agreement with the experimental data when 
(i) the off-diagonal part of the bare Green's function, 
^(k, iw) = Ak/(w^ + E^), is ignored and (ii) the real 
part of A rather then the imaginary part is taken (in 
contrast to the prescription in Eq. (|10|l 'l. The question 
naturally arises: can this prescription be justified in a 
physically well motivated manner? 

(i) may arise simply as a result of chemistry; it is in 
principle possible that impurities in cuprates act both as 
charge and spin scatterers and that the strength in both 
channels is roughly equal. However, such an assumption 
amounts to a fine tuning of parameters. Furthermore, 
chemistry cannot plausibly explain (ii). 

We now show that at least qualitatively both effects 
may arise when the disorder potential is not weak and 
the full T-matrix is considered. In the T-matrix, all of 
the terms in the Born series are resummed by means of 
Eq. Unfortunately, solving the Lippman-Schwinger 
equation for arbitrary disorder potential is a daunting 
task. Progress can be made, however, in the case of 
point-like scaterrers. In particular, for a single point- 
like scatterer at the origin, the T-matrix becomes mo- 
mentum independent and its to dependence can be found 
analyticallyJi One can easily see that the expression for 
5n(q, iv) must be modified as 

5n(q,cj) = -ilmVG°(k,w)f(tj)G'"(k-q,w). (18) 



In the following, we focus on this analytically tractable 
situation which is also relevant to a dilute concentration 
of point scatterers provided that one can neglect coherent 
scattering from multiple impurities. 

First note that in a superconductor, a charge scat- 
terer, which is represented by a potential proportional 
to CT3 in Nambu space, generates terms proportional to 
CTo in the Born series. These terms are magnetic-like since 
the magnetic potential is represented by a unit matrix. 
Magnetic-like terms arise from a pure charge impurity 
due to the 2x2 structure of the Green's function (there 
is no symmetry that forbids their creation). Second, the 
T-matrix is in general complex. Thus, there could be a 
regime in which its imaginary part dominates, implying 
that (5ri(q, to) would receive dominant contributions from 
the real part of A(q, w). 

As shown by Balatsky et al.^ near half filling sum- 
ming the Born series for the point-like charge impurity 
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FIG. 5: The T-matrix approximation for a point like charge 
impurity and a local gap modulation at the origin. The band 
structure parameters and gap magnitude are the same as in 
Fig. □ and y = 0.3t, 5 A = O.SAq. The frequencies from left 
to right are 0.2,0.5 and 0.7Ao. 

yields T{u!) — To(a;)CTo + T3{Lu)a3 with 

W = V ^^(cc') = , ~' . (19) 

where c — cot 5q is the cotangent of the scattering phase 
shift and gol"^) = Tr (k, . The above frequency 

structure of the T-matrix represents a quasi-bound state 
near the impurity with a resonance at c = ±(;o(w). Since 
c is a real number and the Green's function has an imag- 
inary part, at resonance the T-matrix is predominantly 
imaginary and T3(aj) w ^Tq{lj). This means that near 
resonance both (i) and (ii) can be satisfied and the pre- 
scription that we obtained heuristically in the previous 
subsection can be physically motivated. 

Our justification above should be viewed merely as a 
"proof of principle" since the arguments only hold near 
the resonant frequency whereas we need them to apply at 
all frequencies. One could perhaps argue that in real sam- 
ples the impurity resonances will be distributed forming 
a continuous band in which T{lo) is predominantly imag- 
inary and T'i{Lj) « To(tj). However at this stage this is 
merely a speculation. 

Another mechanism that can give rise to suppression 
of the off-diagonal terms is spatial variation of the gap 
amplitude. Such variations have indeed been observed 
in the STS studies on BiSCCO samples.2^ Theoretical 
treatment of these is somewhat complicated by the fact 
that in a dSC the order parameter lives on the nearest 
neighbor bonds of the underlying ionic lattice. For that 
reason we defer this discussion to the Appendix IbI 

Our numerical results for a point-like impurity (in- 
cluding the gap suppression) in the T-matrix approxi- 
mation at the relevant probing frequencies are presented 
in Fig. [SI The obtained results are similar to those ob- 
tained by the simple treatment of section IIC (Fig.^ and 
display most of the features expected by the octet model. 
We emphasize again that considerations in this subsec- 
tion are not meant as detailed modeling of the scattering 
potentials in the cuprates. Rather, they are meant to 
illustrate how the heuristic treatment presented in the 
previous subsection can be reconciled with a more mi- 
croscopic approach. In what follows we shall continue 



to describe the FT-STS patterns 5n(c\,Lo) in terms of 
A(q,tt') via Eq. l(TU|) . 



III. THE PSEUDOGAP STATE 

We now focus on the pseudogap state. In principle, 
the information contained in the FT-STS patterns should 
permit the unambiguous determination of the electron 
order present in underdoped cuprates above Tc (but be- 
low T*). In practice, however, we have seen that details 
of the scattering mechanism can lead to vastly differ- 
ent interference patterns, and even the description of the 
well-understood dSC state is fraught with difficulty. We 
thus set our sights on a more modest goal; in the follow- 
ing we consider the question whether it is possible, using 
the FT-STS techniques, to identify the condensate in the 
pseudogap state as being in the p-p (superconducting) or 
p-h channel. 

As we have seen above the p-p order in the supercon- 
ducting phase gives rise to very special BCS coherence 
factors which mix electron and hole degrees of freedom. 
These coherence factors lead to unique FT-STS patterns 
with sharp peaks at the set of "octet vectors" shown in 
Fig. (Q. According to our analysis there are two cru- 
cial factors which determine the position of these peaks 
in the Brillouin zone. The quasiparticle excitation spec- 
trum Eq defines a set of contours of constant energy, 
Eq = ±u!, along which the denominators of Eq. I|12|) 
vanish and cause large FT-STS response. The coherence 
factors, entering the numerator of Eq. H12|l . determine 
whether the quasiparticle interference is constructive or 
destructive at a given point and thus select special points 
along these contours at which the peaks in FT-STS ap- 
pear. 

It follows from this analysis that the observed FT-STS 
patterns will change when either the excitation spectrum 
or the coherence factors are modified as one passes from 
dSC to the pseudogap state. This observation is the key 
to our proposal to distinguish between the pairing and 
p-h correlations in the pseudogap state. If, on crossing 
the Tc, the FT-STS patterns are essentially unmodified, 
then one can conclude that the pseudogap is dominated 
by ordering in the p-p channel. If, on the other hand, a 
qualitatively different pattern is observed above Tc, then 
this would be evidence for another type of order, most 
likely in the p-h channel. These conclusions should be 
insensitive to the details of the scattering mechanism that 
give rise to the quasiparticle interference. 

A more complicated situation may arise when the com- 
peting order can coexist with superconductivity in some 
region of the phase diagram, as envisioned e.g. in Ref. |^ 
In that case the above analysis still applies but one ex- 
pects a more gradual change of the interference patterns 
as a function of temperature (or another control param- 
eter such as doping) reflecting the gradual growth of the 
competing order at the expense of superconductivity. 

In order to illustrate the crucial role the type of order 
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parameter plays in the process of quasiparticle scattering 
we examine two representative models of the pseudogap. 
Our goal is to understand how the interference patterns 
vary depending on whether the order parameter is in the 
p-p channel (superconductivity) or in the particle-hole 
channel. 



A. Superconducting phase fluctuations 

In this scenario the pseudogap state is viewed as a 
phase-disordered d-wave superconductorii There exist 
several theoretical frameworks to describe such a state 
of electronic matteri2iii2ii2ii4 all sharing a common fea- 
ture that the underlying order is in the pairing channel. 
Here we focus on the "QED3" theorySii^ but it is easy 
to see that other approaches will yield similar general 
outcomes for the FT-STS patterns. 

In the QED3 theory, the phase of the superconducting 
order parameter is disordered by vortex-antivortex exci- 
tations that are encoded in a U(l) gauge field minimally 
coupled to the nodal fermions. At Tc, vortex-antivortex 
pairs unbind through the Kosterlitz-Thouless transition 
and cause a loss of phase coherence at long distances. The 
resulting pseudogap state is described by nodal fermions 
that are strongly interacting, the interaction mediated 
by the massless U(l) gauge field. The fermionic Green's 
function in Nambu spacei^ is given by 

G°(ic...) = A- ,^,;;,+;^-^._„,, . (20, 

where A is a high energy cutoff and rj is the anomalous 
dimension exponent which encodes the dressing of the 
nodal quasiparticles by the gauge field fiuctuations. Gen- 
eral requirements of causality and unitarity dictate that 
77 > 0; however the exact value of the anomalous dimen- 
sion is still under debate>i^iS24^2i2i Most believe that 77 is 
a small positive number. Here we sidestep these issues 
by treating 77 as a free parameter. We show that the 
FT-STS patterns are insensitive to its exact value. 

The similarity between QED3 propagator (|20|l and 
that of the dSC state suggests that the interference pat- 
terns obtained in this model will resemble those in the 
dSC state. We study the FT-STS patterns within the 
same set of approximations adopted for the dSC state, 
i.e. we evaluate A(q, w) and consider its real part. In 
addition, we note that the absence of the off-diagonal 
part in the QED3 propagator (which reflects lack of long 
range superconducting order in the pseudogap state) im- 
plies that the A+A_ term in A(q, w) will automatically 
vanish, irrespective of the type of scattering. Numeri- 
cal evaluations of A(q, w) are carried out using the fast 
Fourier transform technique described above and the re- 
sults are given in Fig. © . As expected based on our gen- 
eral argument the QED3 interference patterns are simi- 
lar to those obtained for the c?SC state. Peaks occur at 
the same points in the Brillouin zone and the anoma- 
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FIG. 6: The interference patterns within the QED3 scenario 
with the anomalous dimension exponent 77 = 0,0.1,0.2,0.3 
and 0.4, ui — 0.45Ao and all other parameters as in Fig. Q. 
The upper panel displays a cut along the (tt, tt) direction of 
the Brillouin zone and the bottom plot is a cut along the (tt, 0) 
direction. In both panels, the solid, dashed-dotted, dashed, 
dotted and solid lines correspond to 77 = 0, 0.4 respectively. 

lous dimension 77 causes some suppression of the sharp 
singularities. 

This similarity can also be understood by perform- 
ing the nodal approximation for the QED3 propagators, 
which is done in appendix ^ This treatment shows 
that the coherence factors produce the same angular de- 
pendence along the elliptic contour as in the dSC state, 
but the square root singularity is replaced by a weaker 
l/(z^ — 1)1/2"'' singularity, where z = uj/\q\ as before. 
Thus, the QED3 and dSG patterns will be qualitatively 
similar as long as < 77 < i, which is exactly the range 
of values where we expect 77 to lie based on general con- 
siderations. 



B. Competing orders 

As an example of a model with competing orders we 
choose the d-density wave (DDW) theoryi& In this model, 
an order parameter of circulating currents with d-wave 
symmetry competes with superconductivity and domi- 
nates the pseudogap phase. The order parameter is in 
the particle-hole channel and connects excitations with 
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FIG. 7: The contour of constant energy for u — 0.3Do in the 
DDW model with /i = ~t, Do = O.lt and t' = -0.3t. The 
dashed Une encloses the reduced Brillouin zone (RBZ) and Q 
is the DDW vector. Note the electron pockets in the vicinity 
of (±7r, 0) and (0, ±7r). 

momentum difference of Q = {tt,7t) (see Fig. ((TJ). The 
Green's function is conveniently written in the basis of 

G°(k, iLj) = [{iLu - e'^) - e'^as - i?k^2]~\ (21) 

with = ^(ek + Ck+q), = iC^k - Ek+q), and the 
DDW gap Dk = i £'0(008 fc^ — cosfcj,). In this model, the 
Born approximation for the density of states modulations 
is given by 

A{q,iLj) = TrJ2'[G°{k,iuj){l + ai)G°{k-q,iu;)] 

k 

{nl + EDini + El) ' ^ ' 

where if2k — ~ ^k ''^^'^ the prime on the summation in- 
dicates restriction to the reduced Brillouin zone depicted 
in Fig. |(7|). Derivation of this simple form is given in 
Appendix 10 It is also worth noting that Eq. H22|l . also 
used by us previously,lSi agrees with that used by Bena 
et alm^ in the appropriate limit of Born scattering. 

When the DDW gap is close to the Fermi energy, i.e., 
close to half filling, the low energy excitations exist only 
near the nodes of the DDW gap. In this case we can 
linearize the spectrum in the vicinity of the nodes and 
approximate: ej^ — > 0, fik ^ (jJ, ^k ~^ vpki and Dk ^ 
VAk2, where ki and ^2 are defined in the same way as in 
the dSC state. The linearized form of Eq. (|22|) turns out 
to be the same as that of a magnetic scatterer in the dSC 
state given by Eq. Hlt)|) . This implies that close to half 
filling, a charge impurity in the DDW model is similar to 
a magnetic impurity in the dSC. This behavior is due to 
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FIG. 8: Interference patterns in a DDW model of the pseu- 
dogap compared with a dSC model. Panels (a-c,g,h) show the 
results in the limit jj. — and t' — 0. Panel (a) represents 
the contour of constant energy at = 0.5Ao. In this limit 
both models have the same CCE. Panels (b) and (c) show 
the calculated LDOS modulations in both models. These FT 
LDOS modulations differ from each other due to the different 
coherence factors as described in the text. These differences 
can also be seen in the one dimensional cut of the Brillouin 
zone along the (7r,7r) direction, shown in panels (g) and (h). 
Panel (d) shows the CCEs for the two different models when 
t' — —0.3t and jj, — — 0.5t. In this case the energy dispersions 
are different; the solid line is the banana shaped CCE of the 
dSC model while the dashed line represents the relevant CCE 
of the DDW model. The calculated FT-STS for this case are 
presented in panels (e) and (f). In all figures ui = 0.5Ao and 
Ao = -Do = 0.5t. 



charge conservation in the DDW model, i.e., there is no 
BCS mixing to produce the angular dependent coherence 
factors of the dSC state. 

The expression H22II can also be evaluated exactly nu- 
merically for all dopings, and we have done so with the 
same band structure as in the previous cases. The results 
are presented in Fig. ((SJ along with the results for a dSC 
with the same parameters. 

The analytical and numerical calculations of the LDOS 
modulations in the DDW model show two principal dif- 
ferences between this theory and the BCS theory of the 
dSC state. First, the very different coherence factors 
lead to arc like patterns rather then peaks in the Bril- 
louin zone. This feature cannot be altered by inclusion of 
higher order terms in the Born expansion (the T-matrix) 
due to the conservation of charge. The solution for the 
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T-matrix equation for the potential V{1 + cti) gives 

fH = (l + ai)- ^^j— (23) 

1 - gQ[uj)V 

where go{(jj) = Tr G°(k, w). This does not affect the 
matrix structure of the scattering vertex but merely adds 
some frequency dependence to it. The second difference 
is in the positions of features in the Brillouin zone. The 
positions are determined by the single particle disper- 
sion, which in the DDW model is similar to that of a 
dSC state only close to half filling and in the absence of 
second nearest neighbor hopping t' . Upon choosing more 
realistic band parameters the CCE of this model develop 
electron pockets and deviate from those of the dSC state, 
as shown in Fig. ^d). This is a result of the fact that 
the node of the DDW order parameter is not attached to 
the Fermi surface but to the ( f , f ) point in the Brillouin 
zone. 

The patterns are most similar in the special case of 
/i = = 0, i.e. exactly at half filling and with near- 
est neighbor hopping only. In this case the excitation 
spectra for DDW and dSC are identical. The coherence 
factors, however, are still different and give rise to qual- 
itative differences in the FT-STS patterns. As shown in 
Fig. ||SJ) away from this unphysical limit the patterns di- 
verge even more and it should be straightforward to dis- 
tinguish between them experimentally. This conclusion 
is independent of whether the superconducting state has 
coexisting DDW order or not. This is due to the fact 
that the FT LDOS spectrum in the mixed state is very 
similar to that of a pure dSG as was shown by Bena et 
alM 

Our conclusions regarding the pure DDW state gener- 
ally agree with those of Ref. S^, based on more detailed 
considerations within the T-matrix approach. In partic- 
ular our original finding,^^^ that the FT-STS patterns in 
DDW and dSC states are qualitatively different, seems 
to emerge independently of details of the band structure 
and the scattering mechanism. We argued above that 
such differences are even more generic and will allow for 
discrimination between p-h and pairing origin of the un- 
derlying electron order in the pseudogap state. 



dering phenomenon is unclear at present and two scenar- 
ios involving a hole Wigner crystal^^ and a pair density 
wave2Si2i have been proposed. 

Abs ence o f the quasiparticle interference patterns in 
Refs. I33l34l implies that the methods described in this 
work are not directly applicable to these experiments. A 
more promissing in this respect is the recent experiment 
of McElroy et al.'^^, which identifies 'zero temperature 
pseudogap' spectra in nominally superconducting sam- 
ples of BiSCCO an low temperatures. Such ZTPG spec- 
tra coexist with the more conventional superconducting 
spectra and form characteristic nanoscale patchwork of 
interconnected regions. In strongly underdoped BiSCCO 
regions exhibiting ZTPG almost completely dominate the 
field of view. The interesting aspect of this data is that at 
low energies these ZTPG regions exhibit FT-STS inter- 
ference patterns that look just like those observed deep 
in the superconducting state. At higher energy, however, 
non-dispersive peaks with periodicity close to 4 lattice 
spacing are observed in these same regions. These obser- 
vations indicate that (i) the low energy quasiparticles in 
the ZTPG regions are of superconducting Nambu-Gorkov 
variety, and (ii) there exists an intimate relationship be- 
tween the superconducting order and the new ordering 
phenomenon these materials. 

An interesting question arises as to why are the disper- 
sive peaks absent in the patterns observed in the pseu- 
dogap state. One possibility is that there are simply 
no sufficiently coherent quasiparticles that could form 
the interference patterns (such as in the QED3 theory 
with 77 > i). Another possibility is that the patterns 
are present but have not been so far detected. It is en- 
tirely possible that in the BiSCCO experiment'^^ done at 
T ~ lOOK delicate interference patterns are washed out 
by thermal broadening. The Na-CCOC experiment)2i on 
the other hand, is done at T ~ O.IK and thermal smear- 
ing cannot be invoked. In this case it is interesting to 
note that no dispersive patterns have been observed in 
this material even in the superconducting state at higher 
doping. These considerations thus raise hopes that the 
interference patterns may be observed in future experi- 
ments and could help determine the nature of the pseu- 
dogap state in cuprates. 



C. Relation to experimental data 

Recently, two groups succeeded in obtaining the 
FT-STS data in the pseudogap state of two different 
cuprate superconductors. Vershinin et al.^- focused on 
the T > Tc region in weakly underdoped samples of 
BizSrzCaCuaOs+a (BiSCCO), while Hanaguri et alM. 
studied lightly doped (non-superconducting) single crys- 
tals of Ca2-2:Na;j;Cu02Cl2 (Na-CCOC) at low temper- 
atures. In both cases no dispersive features could be 
identified. Instead, non-dispersive patterns with period- 
icity close to 4 lattice spacings have been reported in 
both experiments and interpreted as evidence for under- 
lying static charge ordering. Precise nature of this or- 



IV. CONCLUSIONS 

We have shown that the quasiparticle interference pat- 
terns seen in the FT-STS experiments reveal signatures 
of both the quasiparticle dispersion and, through their 
sensitivity to the quasiparticle quantum phases, the na- 
ture of the underlying electronic order in the system. In 
particular, the quasiparticle phase determines the coher- 
ence factors, which then determine the basic characteris- 
tic features of the interference patterns. We have demon- 
strated, by general arguments and detailed calculations 
within two relevant models, that the superconducting or- 
der alone produces patterns consistent with the existing 
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experimental data. In particular, order of another type 
(such as DDW) yields qualitatively different patterns, 
even if the quasiparticle excitation spectrum is identical 
to the dSC. 

The above feature is of key importance since it is in 
general difficult to model the details of FT-STS patterns, 
even in the well understood dSC state. This difficulty 
presumably arises from the complexities of the scatter- 
ing processes which in realistic materials are likely caused 
by several distinct mechanisms and must be modeled as a 
combination of magnetic, non-magnetic and pairing po- 
tentials. In addition the matrix elements for the tunnel- 
ing between the sample and the STS tip, ignored in most 
theoretical treatments, are likely to affect the patterns. 

Despite these complexities, the FT-STS carries a 
wealth of useful information. Sensitivity of the interfer- 
ence patterns to the type of ordering implies that even 
though we might not be able to model their every detail 
in a particular state, we are assured that the patterns 
will change as the material crosses a phase boundary to 
a state characterized by a different electronic order. 

Based on this insight we have proposed a test for the 
nature of the pseudogap phase in cuprates using FT-STS. 
The phase boundary in this case is Tc{x), i.e. the su- 
perconducting critical temperature as a function of hole 
doping X. If the pseudogap is due to fluctuating SC or- 
der, then the FT-STS patterns above Tc should remain 
qualitatively the same as those below Tc- If, on the other 
hand, the pseudogap is due to an order parameter in the 
particle-hole channel, such as SDW, CDW or DDW, the 
patterns above Tc should be qualitatively different, due 
to their different coherence factors. 

Although we focused here on one example of a particle- 
hole order parameter, the arguments presented above are 
general and can be applied to distinguish any such order 
from superconductivity. In addition, DDW order is per- 
haps the most germane to our argument since for spe- 
cially selected parameters the excitation spectrum is the 
same as in the dSC. We have demonstrated that even 
in this special case the FT-STS patterns exhibit qualita- 
tive differences and should be easily distinguishable from 
the dSC. This is a result of the quasiparticle charge non- 
conservation in the superconductor which manifests itself 



in the coherence factors of the FT-STS spectra. Other 
types of p-h order, such as SDW and CDW, will generally 
have very different excitation spectra, trivially implying 
different interference patterns. 

By probing the interference patterns below and above 
Tc(x), future FT-STS experiments should be able to un- 
ambiguously discriminate between the remnants of the 
superconducting order and an order in the particle-hole 
channel, thus settling one of the key puzzles in the high- 
Tc superconductivity. 
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APPENDIX A: LINEARIZED APPROXIMATION 

1. Non- magnetic scattering 

Let us first consider the density of states modulations 
resulting from scattering within the same node in the 
linearized approximation given in Eq. (|13|) . The solution, 
Eq. (|14|) , can be found by "combining the denominators" 
using the Feynman parameterization.^® This is based on 
the formula 



dx-, 



'{1-x) 



A'^BP r(a)r(^) 7o [xA x)B]" 



(Al) 



where r(z) is the Gamma function and the identity is 
valid for any positive A, B, a, (3. Applying this to Eq. H13(l 
with a — P = 1 and performing a shift in the integration 
variable k ^ k -I- (1 — 2x)q yields 



Aiinq,zu; = / dx / — 1 2 , ^2 , a n F2T2 ' 

The ki — fcj term vanishes in the angular integration and the radial integral is convergent. After performing the 
integration over momentum we get 

A / N 1 /■', -^'-4a;(l-a;)(g2_g2) 

Aiin(q,*^) = - dx \ '"-^l^ (A3) 

and the solution quoted in the text follows by performing the x integration and Wick rotating the frequency iuj lo 
to get the retarded quantity. 

Now consider a scattering process between two antipodal nodes, spanned by the vector Q (shown in Fig. Q) 
connecting the nodes along the (tt, tt) direction (the length of Q varies with the doping). The linearized form of A is 
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similar to Eq. H14|l with some signs inverted due to the different directions of vp and wa near the two nodes 

Aim(q+ Q,zw) - 



VFVi, J (2^)2 [c^2 + (k - q)2][c^2 + (k + q)2] 

1 r^^-.-+A.ii-.m-qi)^ ^^^^ 



Attvfva Jo ^'^ + 4a;(l — x)c^ 

where in the last line we have introduced the Feynman parameter x and performed the momentum integral. The 
result is given by Eq. (fTl|l with the angular factor replaced by q\l(f- This 90° rotation of the angular factor causes 
the high intensity regions to be two points on the elliptic contour in the direction parallel to the vector Q as shown in 
Fig. ^p). This is a deviation with respect to the octet model. Instead of the vectors qs and q4 connecting the edges 
of the antipodal CCE bananas we obtain qa « qa at the points Q ± The octet vector q4 is close to the point 

where the angular coherence factor vanishes and hence appears as an end of a line rather than a peak. 

In the two cases above we have scaled the momenta fci and k2 by and wa respectively to achieve radial symmetry. 
This convenient scaling can not be applied to the case of adjacent node scattering. For generic anisotropy one then 
obtains more complicated expressions which cannot be evaluated in a closed form. To keep things simple we treat 
this process in the limit of isotropic velocity vp = v/^ = v. If we denote by Q' the vector that connects two adjacent 
nodes we have 



1 f (fk 

(2^ [c^2 + (k - q)2] [c.,2 + (k + q)2] ' 



Aiin(q + Q ,1^) = — / ,o-^2 r..2 , ;;^2i r, ,2 , , ;;^21 ' (^5) 



where we have omitted terms that are odd functions of 
Lo. These would cancel later, when taking the appro- 
priate components of A according to the prescription in 
Eq. l(Tn|) . The solution is given by 

S(.z) = a.rrtan (A6) 

There is no angular dependence and the resulting A 
is just a uniform contour, as shown in Fig. ^^). 
More detailed considerations indicate that restoring the 
anisotropy would introduce a four-fold modulation of the 
intensity along the contour. 

2. Magnetic scattering 

The inter-nodal magnetic scattering can be calculated 
in a similar way for a magnetic point like impurity, where 
the only difference is the sign in front of the A+A_ term 
in Eq. (|12|l . This leads to the linearized form in Eq. H16|l 
from which it can be readily seen that no angular depen- 
dent coherence factors emerge in this case. The solution, 
Eq. H17|) is obtained by following the same steps as in 
the previous subsection. For antipodal nodes the result 
is the same as for the intra-nodal magnetic scattering, 
Eq. (|12|) . where there is no angular dependence. How- 
ever, when considering the case of scattering between 
adjacent nodes, the magnetic scattering exhibits angu- 



lar dependence that generates peaks close to the octet 
model vectors qi, qs and q2,6 even in the isotropic limit. 
The solution for this case is given by: 

A^ag/ ^ 1 /9l92 + w2\ 

V^q,-) = (A7) 

where S{z) is defined in Eq. (|A6() and the factor of 
9192/9^ is responsible for four peaks along the contour, 
near the octet vectors qi , qs , and q2 e as shown in 
Fig. ©). 

The above analysis suggests that a combination of non- 
magnetic and magnetic scattering is needed to produce 
sharp peaks at all the octet vectors. The exception is 
q4 which marks the end of a line produced by Eq. HA4|I . 
The linearized approximation results are given in Fig. 
In panel (e) we have also combined together the results 
for non-magnetic scattering to cover the entire Brillouin 
zone. For adjacent node scattering we use magnetic re- 
sults which show the four-fold modulation expected for 
the anisotropic case. 



3. QED3 

The nodal approximation is also useful for the case of 
propagators with anomalous dimensions, as in the case 
of QED3 . Let us consider the intra-nodal scattering pro- 
cesses within the QED3 framework. We use the same 
Feynman parameterization as before, where this time 
a = /3 = 1 — 77/2. This leads to the following integral 
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A^(q,iw) = / 

VpVA J 



1 f cPk {iuj + ki + qi){iuj + ki — qi) 



VFVA J (27r)2 [cj2 _^ + q)2]i-V2[,^2 _^ - q)2]i-V2 

™ (r(l - ?7/2))2 "^^J i (2^)2 [c,2 + k2 + 4a;(l-x)q2]2-. 

where the last line was obtained by the variable shift k — > k— (2x— l)q. We can now evaluate the momentum integral 
by first performing the angular integration and then the radial part with a momentum cutoff A. We obtain 

(27r)2 47r(l-77) Vk2 + 4a;(l-a;)q2]i-'' 2??^ \ ^ \ )^ ) \ j \ ) 

1 A. .2 , ^x2lr, 4x(l-a;)gi , 1 



47r(l - r\) 

where only the middle term is divergent. We therefore proceed considering the divergent term only, 

-^'-('^'-)'-"- (i-,)4.....(r(i-,/2))2 1 -^-[-(^-)]"'^ ^.2^4 ^^'"^ 



By performing the Wick rotation and using = uj'^/q^ 
we find 



A^(q,w) cx ( ^ 



da; r , 
[2^(1 



x(l-a;)]i-''/2 



.x)-(f)2]i-r 



(All) 



The angular dependence of A along the elliptic contour is 
already explicit in the above expression and is the same 
as in the dSC intra-nodal scattering case. Unfortunately, 
the Feynman parameter integration can not be performed 
analytically for all values of z. For the case of z < 1 
Eg. lAlll can be evaluated using hypergeometric functions 
but no such expansion exists for z > 1. We therefore 
choose to analyze the integral and find its degree of di- 
vergence. First note that the denominator is divergent 
whenever z'^ = 4x(l — x). For any given value of z'^ this 
leads to two singularities at x/2 = — 1± Vl — z^- How- 
ever, since rj is positive and smaller than 1, the singular- 
ities are integrable unless they overlap. Therefore, the 
integral is only divergent when z2 = 1. In order to find 
the degree of divergence we shift the integration variable 
X ^ X — 1/2 and denote 6 = (1 — z^)/4. The x integral, 
I{z), is given by 



I{z) - 2 



(1 _2:2)i-'//2 



^^"(6-0:2)1-. 

l-r, 



(A12) 



We can now define a new integration variable y = 
and write 

The remaining integral is convergent in the limit & — > 0. 
The integration thus can be performed on [0,oo) and 



the singularity, given by the prefactor, is of the order 
of 1/(1 - z2)i/2-r;^ as stated in the text. 



APPENDIX B: GAP SUPPRESSION IN A 
d-WAVE SUPERCONDUCTOR 

In order to take into account the gap modulations in 
the sample we introduce a local perturbation which will 
be referred to as the gap suppression. We start with a 
lattice model where both the hopping amplitudes t, t' 
and the gap function A^^ are defined on the bonds. We 
assume that A^j = ±Ao on nearest neighbors and is zero 
everywhere else, with plus and minus sign corresponding 
to x and y oriented bonds respectively. On top of this 
uniform state we introduce a gap suppression ^A(R) on 
four bonds around the point R, as shown in Fig. I^. 
The distribution of these impurities is assumed to be a 
featureless function of R. The gap modulation produces 
the following perturbation: 

60. = Yl ['^A(R)44,x(R-R') +h.c.] 

(R,R') 

= ^[<5ApXkC^4-Hh.c.], (Bl) 

where (.., ..) denotes the nearest neighbors on the lattice, 
x(R — R') is the sign of the perturbation, depending 
on the bond orientation and Xk = cosfc^; — cosfcy. To 
leading order, this potential leads to the following LDOS 
modulation 

5n(q,cj) = -i|^Aq|Im[An(q,t^) + A22(q, -tj)] , (B2) 



with 



A(q,cj) 



^G°(k,c^)XkaiGO(k-q,c^). (B3) 

k 
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5A 



(0,0) ^(1,0) 



5 A 



o o o 

(0,-1) 



FIG. 9: The gap modulation local perturbation around the 
lattice point R — (0,0). 



Note that the perturbation Xk'^i is only a function of one 
momentum index, conjugate to the relative coordinate 
R-R'. 

One can go further and find the T-matrix for a com- 
bined local potential of a charge impurity and a gap 
suppression. For simplicity, we assume that the same 
dilute array of impurities is responsible for both per- 
turbations and treat the single impurity at the origin 
in the T-matrix approach. For the combined potential 
— Va^ + SAxi^ai the T-matrix equation © can be 
solved using the ansatz 

T'kk'(t^) = To(cj)cto + T3{uj)a3 + Ti(a;)xkO'i. (B4) 

The frequency dependent coefficients are given by: 



TiH 



V 



= SA 



with 



Aiio) = [1 - gi{Lu)V]^ - [go{Lu)V] 
-gd{uj)5A[l- gi{uj)V]. 



(B5) 



gi{Lo) are functions of the frequency, resulting from the 
momentum integration of different parts of the unper- 
turbed Green's function: 



5o(w) 



k 

E 



Ck 



(B6) 



This T-matrix has a resonant structure in frequency, sim- 
ilar to that found by Balatsky et al. for the case of a 



charge impurity.^- However, we are interested in a partic- 
ular regime of frequencies and scattering strength which 
may not coincide with the resonance. The advantage of 
this formalism in this case is the mixture of different scat- 
tering components and the finite imaginary part of the 
T-matrix. 



APPENDIX C: QUASIPARTICLE SCATTERING 
IN D-DENSITY WAVE MODEL 

In order to calculate A(q, lS) in the DDW model one 
has to carefully define quantities in the reduced Brillouin 
zone (RBZ) shown in Fig. {Tj). Our starting point will be 
the unperturbed DDW Hamiltonian^ 

n = ^[ekfl^flk + iT'kaj^ak+Q -iT'taj^+Qflk] 

k 

^ E'^k^kV'k (ci) 

k 

where ak are the electron operators, the prime denotes a 
sum over RBZ wave vectors and is the DDW order 



parameter. We have defined V'k = (cj;., ^k) such that Ck is 
the usual electron annihilation operator, Ok when k is in 
the reduced Brillouin zone and dk+Q is the usual electron 
annihilation operator when k is outside of the RBZ (and 
k + Q S RBZ). The operators c and d are defined only in 
the RBZ and therefore their labels should be understood 
modulo Q ~ (tt, tt). In what follows we maintain the 
labels of c, d and ip explicitly in the RBZ. In this basis 
the matrix Tik is given by: 



Wk = 



Ck 



fik+Q 



(C2) 



For simplicity we consider point scattering described by 
the Hamiltonian 



6n 



(C3) 



p,p' 



P:P' 

The same conclusion can be derived for arbitrary weak 
potential. 

In order to calculate the tunneling amplitude we use 
the interaction representation, where the DDW Hamil- 
tonian Ti. represents the uniform system and STi. is the 
interaction. We express the Green's function in terms of 
RBZ momenta and expand it to first order in the inter- 
action. The single particle Green's function is given by 



G{r,r':t-t') = -i 



.{Ta^{t)al,{t')Sioo,-oo)) 



{S{oo, —oo)) 



(C4) 



where the S'-matrix is responsible for the time evolution 
of operators by the interaction STi., the angular brack- 
ets denote the vacuum expectation values and T is the 
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time ordering operator. We shall now focus on the equal transform with respect to r 
point propagator G{r, r;t — t') and consider its Fourier 



The momentum indices are defined in the full Brillouin zone, where one of the following four cases may occur. Both 
k and k - q are in the RBZ, both are outside of the RBZ, k e RBZ and k - q ^ RBZ or k ^ RBZ and k - q e RBZ. 
In the RBZ we can write this as 

G(q;t-0 = ' ,, E'(^i[^kW4-q(^0+^kft)4-qft')]e(k-q) + 

[CkW4-q+Q(^') + «fkWcLq+Q(i')]0(k - q)}5((^, -^)) (C6) 

Tr5]'(nV'kWV'Lq(^')0(k- q) +aiVkWV'Lq+Q(i')0(k- q)}(i')^(oo,-(X3)). 



—I 



(S'(cx),-oo)) ^ 



where 



QW = {J: MRBZ' 0(k)-i-e(k). 

In the last line of HC6|I we introduced the ai Pauli matrix in the trace in order to implement the sum over off-diagonal 
elements of the 2x2 matrix V'kV'k-q- 

We now include the perturbation in the S'-matrix to first order, i.e., we approximate S'(oo, — oo) w 1 — i J STidt. 
Dropping the unperturbed part of the Fourier transformed Green's function we find the first order correction 

G(i)(q; t-t')^ -VTtJ2' J dt"{T{Mt)4'L^it')e{\^ - q) + ^iVk WV'Lq+Q(i')0(k - q)}V'^r)(l + ^i) 

kpp' 

(C7) 

Performing Wick contractions and a time Fourier transform yields the quantity A(q, cj), defined in the text: 
A(q,w) = ^ J dte"'^'^'-*'^G^^\q;t-t') 

= Tr^' f e-"^(*-*')[G"(k,i-t")(l + (Ti)G*'(k-q,t"-<')0(k-q) 

k 

+ criG°(k, t - t")(l + cri)G°(k - q + Q, i" - i')0(k - q)] (C8) 
= Tr^'[G"(k, w)(l + CTi)G°(k - q, uj)e{k - q) + aiG"(k, + CTi)G"(k - q + Q, w)e(k - q)] 

k 

where G°(k, w) is the unperturbed DDW Green's function given in Eq. (|21|) . The above expression can be simplified 
further. Let us concentrate on the second term when k — q ^ RBZ. If wc now extend the definition of the Green's 
function to hold even when the argument is not in the RBZ (using the functional definitions of ek and -Dk), we find 
that G"(k + Q,Ci;) = (TiG°(k, w)(Ti. Rewriting the second term we find: 

Tr[c7iG°(k, Lu){l + ai)G°(k + q + Q, c^)] = Tr[aiG°(k, lu){1 + ai)aiG"(k + q, Lu)ai] 

= Tr[G°(k,c^)(l + ai)G"(k + q,w)]. (C9) 

Therefore 

A{q,u;) - ^'Tr[G"(k, w)(l + ai)G°(k + q,iv)], (CIO) 

k 

where k + q need not be in the RBZ. 
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